Bistability and noise-enhanced velocity of rolling motion 
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We investigate the motion of a hard cylinder rolling down a soft inclined plane. The cylinder 
is subjected to a viscous drag force and stochastic fluctuations due to the surrounding medium. 
In a wide range of parameters we observe bistability of the rolling velocity. In dependence on 
the parameters, increasing noise level may lead to increasing or decreasing average velocity of the 
cylinder. The approximative analytical theory agrees with numerical results. 
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PACS numbers: 46.30.Pa, 62.40.+i, 81.40.Pq 

Since 1785 when Vince described systematic experi- 
ments to determine the nature of friction laws the 
effect of rolling friction has been investigated by many 
scientists according to its great importance in engineer- 
ing and natural sciences (e.g. [2( and many references 
therein). 

In this letter we study the rolling motion of a hard 
cylinder on an inclined soft plane. Experimentally it was 
observed that the rolling friction coefficient depends non- 
monotonously on the velocity Q: For small velocities 
the rolling friction force increases with increasing veloc- 
ity, while for fast motion it decays as the velocity grows. 
Furthermore, the cylinder is subjected to the drag force 
of the surrounding medium, e.g. air, and unavoidably 
present fluctuations. We will show that this system re- 
veals non-trivial dynamics, i.e., bistability of the cylin- 
der's velocity and noise controlled average rolling veloc- 
ity. 

Model. We consider a cylinder of radius R, mass M 
and moment of inertia I which rolls at velocity v on a 
plane Q. The cylinder is subjected to an external driv- 
ing force F cx — Mg sin a due to the plane's inclination by 
the angle a. The rolling friction force Fr and the viscous 
drag force Fr> due to the surrounding air counteract this 
motion. We also assume that the tangential force act- 
ing between the cylinder and the surface at the contact 
area is strong enough to keep it from sliding. For certain 
materials it has been shown that surface effects such as 
adhesion may have significant influence on rolling friction 
(e.g. 0). On the other hand, for viscoelastic materials 
it was reported that rolling friction is due very little to 
surface interactions, i.e., the major part is due to defor- 
mation losses within the bulk of the material Q. Here 
we assume that surface effects may be neglected. 

Hence, Newton's equation for the cylinder reads 

(M + I/R) v = -F D {v) - F R (v) + F cx + ({t) . (1) 

The stochastic force £(t) describing fluctuations in the 
media is modeled by Gaussian white noise of zero aver- 
age, (£(*)) = 0, and intensity a: (£(*)£(*')) = a 2 S(t-t'). 

For the rolling friction force Fr we adopt the model 
developed in Q, where it is assumed that the plane is 




FIG. 1: A rigid cylinder of radius R rolls down an inclined 
plane. The cylinder contacts the surface of the deformed plane 
in the region (£_,£+) which moves along with the cylinder. 



composed of springs of elasticity k, mass m and damping- 
constant 7 to model elastic, inertial and viscous prop- 
erties of the surface material. The rolling friction orig- 
inates due to the dissipative and inertial resistance of 
the springs, which are compressed and subsequently de- 
compressed when the cylinder rolls on the surface. As 
an important property this simple model leads to a non- 
monotonous rolling friction as a function of the velocity 
as observed in experiments 0. Let us briefly sketch 
the main results of the model which will be used below. 
The rolling friction force is given by 
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where £ is the coordinate in the comoving frame (Fig. [T]l. 
/(£) is the corresponding force density exerted from the 
compressed material on the cylinder, 
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and h(v) is the deformation. The force /(£) acts in the 
contact contact area < £ < £,+ , where the front bound- 
ary 



£+ = V2Rh 



(4) 



follows from geometry, while for the rear boundary the 
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condition /(£_) = yields 



£_ = 7u/fc - y / 2hR+{-/ 2 /k 2 ~2m/k)v 2 



(5) 



The deformation may be found from the balance of forces 
normal to the plane, 
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The first term in Eq. (J2J) equals, up to the factor 1/R, the 
total torque exerted on the body due to its interaction 
with the deformed surface (£/(£) is the torque density), 
while the second term in Eq. @ as well as in Eq. © 
accounts for the surface-cylinder interaction at the in- 
finitesimal region at £ + where the acceleration diverges 
as the cylinder rolls over (for a detailed discussion of this 
peculiarity see Q). In this approach it is also assumed 
that the rolling body does not induce surface waves. 

The viscous drag force due to the surrounding air is 
given by 



Fr 



-Av - Bv 2 



(7) 



where the (positive) coefficients A and B depend on the 
Reynolds number Re, which characterizes the motion of 
the body. These coefficients may be estimated from hy- 
drodynamic theory ||. 

For cylinder velocity ~ 1 — lOm/s, radius ~ 0.1 — 1 m 
and kinematic viscosity of the air v — 0.15 TO -4 m 2 /s the 
Reynolds number is Re ~ 10 4 — 10 5 , which corresponds 
to the drag force Fjj — 1 \cpv 2 S with the constant c ss 0.5 
(here p is the air density, S - is the body cross-section) [ij . 
This yields the estimate for the constant B ~ (0.01 — 1) 
kg/m 2 . 

For very small velocities, i.e. for Re <C 1, the 
drag force is approximately proportional to the veloc- 
ity. Lamb's relation for a cylinder (Fn = Airvpv/b with 
b = log [3.70 z;/ (Re v)], Rev/v < 1 H ) yields the esti- 
mate for the linear constant, A < 2 • 10~ 4 kg/s/m. Thus, 
for the parameters addressed here we can neglect the lin- 
ear drag force, i.e. A = 0. We studied the above model 
analytically and numerically and observed the overmen- 
tioned effects of bistability and noise controlled velocity 
for a wide range of parameters. Here we report the re- 
sults for a cylinder of mass per unit length M = 100 
kg/m and radius R = 0.5 m. The surface material pa- 
rameters correspond to soft rubber 0: m = 100 kg/ m , 
k = 10 6 - 10 7 kg/m 2 /s 2 , 7 = 10 4 '- 10 5 kg/m 2 /s (all 
parameters per unit length). 

Analytical approximation. As shown below for the 
given set of parameters the effects of interest, namely 
bistability of the rolling velocity and noise controlled ve- 
locity, occur at rather large values of the velocity. There- 
fore, as a simplification we assume that the material does 
not recover immediately after the cylinder has passed by, 
but with a delay, when the cylinder has already left. The 
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FIG. 2: a) Rolling friction force Fr and b) deformation h(v) 
over the velocity v. The full lines show the numerical solu- 
tion of Eqs. 1415161 while the dashed lines show the approx- 
imative theory, Eqs. JHJ and @, respectively. Parameters: 
m =100 kg/m 2 , k = 10 7 kg/m 2 /s 2 , 7 = 5 x 10 5 kg/m 2 /s, 
B = 0.2 kg/m 2 . 
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FIG. 3: Total force F(v) for different angle of inclination (the 
number attached to the curves is the inclination in degrees). 
Parameters: (a) as given in Fig. [5] (b) same, except for 7 = 
0.8 x 10 5 kg/m 2 /s, k = 10 6 kg/m 2 /s 2 , B = 0.65 kg/m 2 . Full 
lines: numerical solution, dashed lines: approximative theory. 



according approximation £_ = means that the rear part 
of the material does not recover on the time scale of the 
cylinder motion ~ R/v. A similar approximation was in- 
vestigated in . The balance equation JBJ reduces then 
to Mgcosev = jvh + 3 \/2i?fc/i 3 / 2 , which yields with the 



condition of large damping, vj 3> ky / 2Rh, the velocity 
dependent deformation 



h(v) 



Mg cos a 



1i / 2RMg cos a 
3 V 7« 
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With the same arguments from Eq. we obtain the 
rolling friction force 



ktf 



mv h 
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Bistability. Solving Eqs. 141516(1 we obtain the de- 
formation h(v) as a decaying function of the velocity v 
(Fig. [2h), he., the faster the cylinder moves, the less 
the surface is compressed. Using the obtained deforma- 
tion h(v), the corresponding rolling friction force may be 
calculated (Fig. Eli). Fr depends non- monotonously on 
the velocity, first it increases while at larger velocity it 
decreases. The non- monotonous behavior of Fr(v) re- 
sults from the interplay between the amount of deformed 
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FIG. 4: Bifurcation diagram for stable velocity states (full 
line: numerical solution, dashed line: analytical approxima- 
tion), (a) corresponds to Fig. [5] and Fig. [3^, (b) corresponds 
to Fig. [3Jj. In both cases the agreement between numerics 
and theory is within few percent. 




FIG. 5: Force potential U(v) (F(v) — —dU/dv) over the ve- 
locity v (parameters as given in Fig. |5J for angle of inclination 
1°, 2° and 3°. 



material and the deformation rate: For small velocity, 
Fr(v) increases since the surface material is compressed 
at increasing rate with growing v while at higher rolling 
velocity the amount of deformed material decreases as 
h(v) decays. Consequently for fast motion the rolling 
friction force decays with increasing velocity. 

Given the non-monotonous rolling friction force, the 
total force acting on the cylinder, T(v) = Fo + Fr + F cx , 
is represented for realistic parameters by a iS-shaped 
curve, implying bistability of the velocity for a certain 
range of parameters. Indeed, the steady state condition, 
v = 0, or T{v) — may be fulfilled cither for only one 
velocity (top and bottom curves in Figs. [3]or for three dif- 
ferent velocities (middle curves). The former case implies 
a unique stationary velocity, while the later case allows 
for three velocities. Only two of them, the smallest and 
the largest, correspond to stable motion (dT/dv < 0). 
By varying the parameter of the drag force B, or by mod- 
ifying the driving force M g sin a the system may transit 
from one stable solution to the other. Figure 0] shows the 
bifurcation diagram. 

Starting from an inclination angle that corresponds to 
a monostable velocity, the variation of a causes the ap- 
pearance of the bistable regime, change of the symmetry 
of the potential U(v) (Fig. |3J and the transition to the 
other monostable regime. 

As shown in Figs. 121 II the analytical theory (dashed 
curves) agrees well with the numerical results if the 
damping parameter 7 is large. The deviation in Fig. 0] 
appears large due to the very detailed resolution of the 
horizontal axes, in relative units the deviation is few per- 
cent. 

Noise-induced jumps. Up to now we did not consider 
noise, which is always present in realistic systems. If 
the system has only one stable velocity, the addition of 
noise does not change the motion of the cylinder qualita- 
tively. In this case the velocity fluctuates around the av- 
erage value given by the steady state condition T{v) = 0. 
Figure HO (top left and right) shows the velocity of the 
cylinder as obtained from a numerical integration of the 
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FIG. 6: Top: velocity of the cylinder subjected to Gaussian 
white noise of intensity <r 2 =50 kg 2 /s 3 as obtained from the 
numerical integration of Eq. ijljl. Bottom: the corresponding 
velocity distributions of the cylinder. From left to right: a = 
4.87°, 4.9°, 4.93°, all other parameters are the same as in Fig. 
3b. 

stochastic equation Q). The corresponding velocity dis- 
tribution reveals a single peak (Fig. bottom left and 
right). 

For parameters which correspond to bistable velocity 
the presence of the noise changes the system behavior 
qualitatively. Instead of moving with a fixed velocity 
there are stochastic jumps between the meta-stable ve- 
locities. Figure HJI (top, middle) illustrates this regime. 
Consequently, the velocity distribution reveals two well 
separated peaks (Fig. bottom, middle). 

Noise-enhanced velocity. Under certain conditions due 
to the nonlinear dependence of the rolling friction on 
the velocity the average cylinder velocity may increase 
with increasing noise level. Since v increases at the same 
driving force F cx , one effectively has increased mobility 
fi = v/F. This surprising effect is shown in Fig. \7\ 

This phenomenon may be understood analyzing the 
dynamics of the system with respect to the force potential 
U(v). In the absence of noise the velocity remains in its 
(meta-)stable state, where dU/dv = —T(v) — 0, i.e., at 
the minimum of the potential U(v) (Fig. 0). Subjected to 
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FIG. 7: Mean velocity of the cylinder over the noise intensity, 
(a) a = 4.87°, (b) a = 4.93°. The parameters of the system 
are given in Fig. |3j5. 

the noise, the velocity of the cylinder fluctuates around 
this minimum. If the potential U (v) has two minima, 
the velocity jumps between both potential wells. If one 
potential well of U(v) is narrow while the other is wide, 
the average velocity will be shifted towards the wider 
well. This shift increases with increasing transition rate 
between the wells. If the potential well is wider for the 
larger stable velocity than for the smaller one, the average 
velocity, thus, increases with increasing noise level (Fig. 

In the opposite case, large noise impedes rolling 
(Fig. 03) . 

Summary. We investigated the rolling motion of a hard 
cylinder on an viscoelastic plane in the presence of a sur- 
rounding medium. The motion of the cylinder is driven 
by an external force due to the inclination of the plane. 
The rolling friction force and the viscous drag force due 
to the surrounding medium counteract its motion. For 
certain realistic parameters the stationary velocity of the 
cylinder is bistable. The numerical results agree well with 
the predictions of a simplified analytical theory in the 
limit of large damping. 

In the presence of noise as it is unavoidable in any re- 
alistic experiment, noise- induced transitions between the 
meta-stable velocities are observed. Depending on the 
system parameters, increasing noise level may accelerate 
or decelerate the rolling motion. 

The described effects may be important for technical 
systems where the addition of noise may lead to an effec- 
tive increase of the mobility of a rolling body driven by 
an external force. 
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